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Abstract 

A medium modified gluon propagator is used to evaluate the scattering cross 
section for the process gg — ► gg in the QCD medium by performing an ex- 
plicit sum over the polarizations of the gluons. We incorporate a magnetic 
sreening mass from a non - perturbative study. It is shown that the medium 
modified cross section is finite, divergence free, and is independent of any 
ad-hoc momentum transfer cut-off parameters. The medium modified finite 
cross sections are necessary for a realistic investigation of the production and 
equilibration of the minijet plasma expected at RHIC and LHC. 



PACS: 12.38.Mh; 14.70.Dj; 12.38.Bx; ll.10.Wx 



Typeset using REVTeX 



*present address: T-8, Theoretical Division, Los Alamos National Laboratory, Los Alamos, NM 
87545, USA 

1 



A lot of research is undertaken to detect a new state of matter, known as the quark-gluon 
plasma (QGP) which is believed to have existed in the early stage of the universe, ~ 
seconds after the big bang. |I],|2| Lattice QCD calculations indicate that this deconfmed state 
of quark gluon matter exists at high temperatures (~ 200 MeV) or high energy densities 
(~ 2 GeV j fm 3 ) ||. The relativistic heavy-ion colliders RHIC (Au-Au collisions at a/s = 
200 A GeV) and LHC (Pb-Pb collisions at ^Js = 5.5 A TeV ) offer the unique opportunity to 
study the production of this state of matter in the laboratory jlj. Perturbative QCD (pQCD) 
estimates that the energy density of jets and minijets produced in these collisions could be 
larger than about 50 and 1000 GeV/ fm 3 at RHIC and LHC respectively 0,0. However, even 
if the attained energy density is sufficient to form a quark-gluon plasma, it is not at all clear 
whether the parton medium will equilibrate and form a thermalized quark-gluon plasma 
before hadronization. Evidence for the formation of a QGP at RHIC and LHC can only be 
established if various proposed signatures have been examined experimentally. The most 
prominent signatures suggested so far are J/^ suppression H], strangeness enhancement 



0, and dilepton and direct photon production |pHTl|l. A quantitative calculation of these 
predicted signatures will be difficult to achieve for the situations found at RHIC and LHC. 
One main uncertainty originates from the absence of a reliable study of the space-time 
evolution of the parton distributions in these experiments. An investigation of this evolution 
will determine the equilibration times and the time evolution of all other local and global 
quantities such as energy densities, number densities, and temperatures of the SU(4)-flavour 
components of the quark-gluon plasma. 

The space-time evolution of the QGP during an ultra-relativistic heavy-ion collision 
might proceed through different stages such as: 1) the pre-equilibrium, 2) the equilibrium, 
and 3) the hadronization stage. During the equilibrium stage, hydrodynamics can be used to 
describe the dynamics of the QGP. However, both the pre-equilibrium and the hadronization 
stages have to be studied in more detail in order to see whether either partons or hadrons 
equilibrate or not. The first stage starts just after the initial nuclear collisions, where many 
hard, semihard, and soft partons are produced. The hard and semihard partons (jets and 
minijets) formed at RHIC and LHC can be desribed by using pQCD [fHf)Ul2| ) , and soft 
gluons may be described by the creation of a coherent chromofield [|I3|-|i7||. For the sake 
of simplicity, we restrict our investigation to jets and minijets only. The equilibration of a 
minijet plasma can be studied by solving the relativistic transport equations with binary 



and secondary parton-parton collisions taken into account [II 

p»dj(x,p) = C(x,p). (1) 

In equation ([[]) 

° { - X ' P) = I {Lfp% {Lfpl {Lfp\ l M (^ 2 ^P3P 4 )| 2 [f(x, P z)f(x,p 4 ) - f(x,p)f(x,p 2 )} 

(2) 

is the collision term for the 2 body partonic scattering process pp<i — > PzP^- The collisions 
among the partons determine how the QGP reaches equilibrium. So, the investigation of the 
2-body collision cross section in the non-equilibrated medium is of great importance for the 
understanding of the minijet plasma evolution at RHIC and LHC. The different partonic 
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scattering processes which have to be considered are e.g. gg — > gg, qq — > gg, gg — >■ gg etc. 
and in more detailed studies also higher order processes like — > ggg might have to be 
taken into account. Most probably, the dominant part of the minijets consists of gluons. 
Hence we consider collisions of the type: gg — > gg. The differential cross section for this 
process is given by: 

da 97va 2 ut us st 

dt = ~2?~ [ ^~ ~¥~rf 1, { } 

where s = (q x + q 2 ) 2 = (q 3 + qif, t = {qi~ q±) 2 = (<?2 - g3) 2 , and u = {qx~ g 3 ) 2 = (g 2 - g 4 ) 2 
are the Mandelstam variables of the partons. For real gluons, they are related by 

s s 
t = - - [1 - cos cm ] and u = - - [1 + cos 9 cm ), (4) 

where 8 cm is the center of mass scattering angle, which goes from — > ^ for identical particles 
in the final state. If \M(s,u,t)\ 2 is put from Eq.(D) into Eq.(0), a divergence emerges in 
the collision term at small angle (8 cm — ► 0), which corresponds to collisions with small 
momentum transfer (t — > 0). Hence, to obtain a finite collision term, one is forced to put a 
minimum momentum transfer by hand, as it is done in the parton cascade model []BJ. The 
plasma evolution and the determination of all the signatures then crucially depend on the 
choice of this momentum cut-off parameter. 

The equilibration time and equilibration process can be studied by using a collision term 
in the relaxation-time approximation for the transport equation |]20|j2T| ,|4]| : 

C(x,p) = [f(x,p) -/ e9 (x,p)]/r c (r), (5) 

where r c (r) is the time dependent relaxation time and f eq = (e pMu,1//T ( T ) — is the 

Bose-Einstein distribution function, represents the flow velocity. The time dependent 
relaxation time t c (t) can be determined from the transport cross section (a tr ) and the 
number density (n) of the plasma via the relation jl[20^2|: 



Tc(r) = — 1 . (6) 
a tr (T) n(T) 

Here, the time dependent number density is determined from the non-equilibrium distribu- 
tion function via: 



(7) 



n(r) = 

and the transport cross section for the gg — > gg collisions is given by: 

a t (r) = £,dt d ^ sm 2 6 cm = f dt ^ \M(s,t,u)\ 2 ^. (8) 

Hence, whether one uses Eq.(H) or (|^) as the collision term, one always encounters the 
divergence in the limit t — > 0. 

This divergence is inescapable as long as we are considering the propagator in the vac- 
uum. As the equilibration of the quark-gluon plasma crucially depends on the cut-off used 
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to remove these divergences JT9|, it is essential to incorporate the medium effects to ob- 
tain finite and cut-off independent partonic scattering cross-sections in the QCD medium. 
Here, one should therefore use the the medium modified propagator instead of the vacuum 
propagator. In this paper, we evaluate and analyze the matrix element squared, total cross 
section, and transport cross section for the process gg — > gg by using the medium modified 
propagator performing an explicit sum over the polarizations of the physical gluons. This 
analysis is important to study the evolution of the QGP because the divergence is removed 
automatically and cut-off independently. 

In this paper we choose the medium modified gluon propagator in the covariant gauge. 
It can be split into the longitudinal and transverse parts [2^,24]]: 

N*"{K) = ^ + ^ + (a - 1)^-4- (9) 

V ; —K 2 + n T -K 2 + U L V ' K A y ' 

where a is the gauge fixing parameter. P£ v and Pj* u are longitudinal and transverse tensors 
given by: 

_ -(uK* - K 2 U> 1 ){uoK v - K 2 U V ) 

Pj 



L 



K 2 k 2 



Pf = gT + ^-, (io) 

9fiu = 9fiu — UyU v with Up being the flow velocity. The quantity u> = K • U is interpreted 



as the Lorentz- invariant energy and k = y—K^K^ with = — U^(K ■ U) as the 

three momentum of the virtual boson. In the local rest frame (U = (1, 0)), cu and k are its 
energy and momentum. The expressions for IT^ and Ut in the high temperature expansion 
are derived by perturbative methods in [23]: 

Ilr = m 2 n (l — x 2 )\l log I X I + i—x] 

n T = mli^ + |(1 - x 2 ) log |1±|| - ijx(l - x% (11) 

where x — x and mo = g 2 T 2 . The t-channel matrix element for the process gg — > gg in the 
medium is given by: 

M = g 2 fae d fe b AelV x ™(- qil q x - g 4 , ?4 )A tV (?i - qjV^fa - q 3 , -q 2 , q 3 )e^ 3 , (12) 

where V^ Xu {j>i,p2,Pz) = [(pi — P2) u 9 flX + (j>2 — PzYg Xv + (P3 — Pi) x 9 tJ,u ] is the three gluon 
vertex, e^(p) are the polarization vectors of the gluons and A Ml/ (/c) is the medium modified 
gluon propagator given by Eq. (^|). 

To sum over initial and final spins of the gluons in the matrix element squared, we use the 



appropriate projection operators for the transverse polarization states of the gluons [25,26 



V- C A *A' _ AA' , 2(glV + 9l 92 ) V- u W _ W , 2 ^^' + 

e i e i — ~9 H 7 — — — ^ > e 2 e 2 — — 9 + 



spins w x 1 spins 



(qi + Q2) 2 ' SDins (gi + aO 1 



E #T = -<T + 7 4 ;7 ; , and £ = -g™' . (13) 



spins ' '/ spins 
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After summing over the final and averaging over the initial spins and color, we obtain (in 
Feynman gauge): 

II 2 A B 

I I medium 



S \t - (t v • uynu L - t)(u L - 1) s*(t - (t v • uff(u T - t)(u L - 1) 

D + ^ n = ^ (14) 



s %t- (t v - umu L -t)(n T -t) s 3(t_(t v .[/)2)2(n T _ t) (n r _t) 



where A, 5, C are: 



A = |t 2 (7 4 {(s 2 [(s v • C/ - u v ■ U) 2 + 4(t - (t v • f/) 2 )] + 2*s [2t + (« v • U) 2 - (s v • U) 2 - 2(t v ■ U) 2 ] + 

o 

2t 2 [{u v ■ U) 2 + 2{u v ■ U)(s v ■ U) - (s v • U) 2 }) x (t [2t + [u v ■ U) 2 - 2(t v ■ U) 2 - (s v • U) 2 } 
+ s [(s v • U + u v ■ f/) 2 + 3(t - (t v ■ U) 2 )])} (15) 



B = -tg 4 {2s 4 [t - (t v ■ t/) 2 ][(n v ■ f/) 2 - (s v ■ Uf] + s 3 t [(s v • Uf + 2(t - (t v • U) 2 - (u v • Uf) 

o 

(s v • t/) 2 + 2{t - (t v • U) 2 ){u v ■ U){s v ■ U) + (uy ■ Uf(U + (uy ■ U) 2 - 4(t v • U) 2 )} 
+s 2 t 2 [-3{s v ■ Uf - 2{u v ■ U){s v ■ Uf + 2(t - (t v ■ U) 2 ){s v ■ U) 2 + (u v ■ U){-5{t v ■ U) 2 
+2(u v ■ U) 2 + 5t)(s v ■ U) + (uy ■ U) 2 (-5(t v ■ Uf + 3K ■ U) 2 + 5t)] + 2st\u v ■ U) 
[2{u w ■ U){s v ■ Uf + 4(-(t v • U) 2 + (u v ■ Uf + t)(s v • 17) + K • C0(-3(* v • U) 2 + 
2(u v • f/) 2 + 3t)] + 2t 4 [s v • [/ + ii v • U] x [(s v ■ £/) 3 - 3u v ■ U{s v ■ U) 2 
+ (u v ■ Uf(s v ■ U) + (uy ■ U)(t + (u v ■ U) 2 - (t v • U) 2 )}} (16) 

C = lg 4 {4s 5 [t - (t v • U) 2 } 2 + 4s 4 t [t - (t v ■ f/) 2 ] [t - (s v ■ f/) 2 - (U ■ U) 2 + K ■ U) 2 } 

o 

+ s 3 t 2 [(s v • Uf + (3(t v • U) 2 - 2( Mv • U) 2 - 3t)(s v • t/) 2 + 2((t - (t v • t/) 2 )(w v • C/)(s v • ?7) 
(uy • f/) 4 + 4(t - (t v • t/) 2 ) 2 + (t - (t v • £/) 2 )K • t/) 2 )] + s 2 t 3 [-3(s v • f/) 4 
-2{u v ■ U){s v ■ Uf + 8(t - (t v ■ ?7) 2 )(sv • f/) 2 + 2u v ■ C/(-3(t v ■ f/) 2 + {u v ■ Uf + 3t)s v ■ 
+3K ■ f/) 4 + 7(t - (t v ■ f/) 2 ) 2 - 6(t - t v ■ U 2 ){u v ■ Uf)] + 2st 4 [2(uy ■ Uf + 4(s v • U)(u v ■ Uf 
+(2(s v • f/) 2 + (t v • Uf - t)(w v • Uf + 2s v ■ U(t - (t v • Uf )u v ■ U + (t - (t v • t/) 2 )(7(s v • ^) 2 
-(t v ■ t/) 2 + t)] + 2t 5 [(s v • f/) 4 - 2u v ■ C/(s v ■ f/) 3 + 2(-(t v ■ Uf - (« v • t/) 2 + t)(s v • Uf 
+2u v ■ U((t v ■ Uf + (u v ■ Uf - t)s v -U + {u v - Uf - (t - (t v ■ Uf) 2 }}. (17) 

U is the flow velocity and s v , t v , w v are defined as follows 

< = (9i + 92)^ = (93 + g 4 r, ^ = (g 3 - = (9i - 94) M , < = (9i - 9 3 r = (94 - 9 2 ) M . (IS) 
The Mandelstam variables of the partons are just given as: 

Eq. ([14]) together with Eqs. (|i~5|) to (|I7|) show the general expression for the t-channel 
matrix element squared of the scattering process gg — > gg in the medium. As this is a 
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rather lengthy expression, we analyze the behavior of this scattering term in the local rest 
frame of the fluid U = (1,0). In the center of mass frame of the parton-parton collision Eq. 
(|H|) becomes: 

,da. _ 97m 2 r (s 2 + 2ts + 2t 2 ) 2 (s + t)(s 2 - 2t 2 ) 

) medium n l Af~n j.\ ffj +'\ 



K dt jmcamm 8 l s*(Tl L -t)(Tl L -t) s 3 (n T - t)(U L - 1) 

(s + t)(s 2 -2t 2 ) (s + t)(s 4 - 3ts 3 + 15tV + 8t 3 s - 2t 4 ) 
+ s ^(U L -t)(Tl T -t) + s 5 (U T -t)(Tl T -t) } - ( } 

Using the in-vacuum gluon propagator and repeating the above procedure one finds for the 
t-channel process: 

da 9ira 2 4s 5 + 4ts 4 + 16tV + 27t 3 s 2 + 10t 4 s - 2t b 

(^)vacuum - — ^ • (21) 

It can be verified that for = and Ht = Eq . (|20|) reduces to Eq . (pT|) . In the limit t — > 
the dominant divergent part of Eq. ( ^Tj ) is given by: 

da _ W 

dt ' eac ^ vacuum ' \ J 

This divergence can be removed by introducing a screening mass mjj by hand pTPfl: 

,da . 9na 2 



dt jlcad.cut-off 2 (t-ml,) 2 " 



(23) 



The introduction of the Debye screening mass (the longitudinal part of the self energy) in 
the above formula leads to a screening of the long range electric field but not of the magnetic 
field. However, Eq. fl2~0|) can be used to screen both the electric and magnetic (beyond one 
loop in self energy evaluation) part simultaneously. For a comparison with Eq. (^) we 
extract the terms quadratic in t in the denominator of Eq. (|20|): 

da 9na 2 1 1 

(^)lead.medium = — { ^ _ f) ^ _ f) + ^ _ f) ^ _ f) 

+ (u L -t)(ii T -t) + (u L -t)(n T -t) } - (24) 

This equation gives the principal contribution to the total cross section (see below). Hence, 
this expression can be used for all practical purposes during the evolution of the plasma. 
In the situation we are considering, it can be easily checked that in the leading order = 
m 2 D = g 2 T 2 and Hp = 0. This suggests that even if we consider Eq.(^4j), the magnetic 
part is still not screened as long as we use the perturbative expression of the self-energies 
in the leading order. In this study, the logarithmic singularity in the transverse part of the 
total cross section is removed by introducing a non-perturbative magnetic screening mass 
m ma g = §(0.255g 2 T) 2 taken from p7jj28|. With these self energies Eq.(f2~4"D can be integrated 
to obtain the total cross section: 



(fltot) lead. medium / 



° A ( — \ 

dt ' ea< ^ mecmum 



W g , 1 , 2 ^Tfc + 2)-ln(^ + 2) 2 

8 y u L ( s + 2U L ) n T n T -n L s + 2n T ; ' 1 ; 
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and the transport cross section: 

r° da 4ut 

\ytx) load, medium J Q>v \ J lead. medium ^2 

97ra 2 r j n (^_ + 2 i) + i n (_^_ + i) 4 l l 



■{2- 



+ 



4 s s s + 2U L s + 2U T 

81n(^r + l)nj - 2( 3 (1 - 21n(^ r + 1)) + 2(ln(^ + 1) - ln(^ + l))n r )II L 

s 2 (n L -n T ) 

2n T (s(l - 2 ln(^ + 1)) - 4 ln(^ + 1)H T ) 



+ _^ ,^; nT) - — >■ (26) 

To show that the above leading order (in t) cross sections give the dominant contribution, 
we evaluate the total and transport cross sections by using the t-channel expression for % 



as given by Eq. (|20|) and then plot the differences between the full and the leading order 
expressions for the total cross section in Fig.l. The graphs show that the difference between 
the full cross sections and the leading order cross sections are very small. So, for practical 



purposes, one can use Eq. (|25|) and (|26|) as the total and transport cross sections for the 
process gg — > gg in the medium. In the figures we use T — and a = 0.3. 

The importance of the medium modified total and transport cross sections given in Eq. 



and fl26|) will be discussed below. To demonstrate it, we evaluate the total and transport 



cross section by using Eq. (p3|) which yields: 

f° da. 9na 2 s 

(Ctot)lead.cut-ofT = dt { — )lead.cut-off = ("T^rt — , o tt ) ' (27) 

J-| at + ZsLi-L 



and 



( Ctr) lead, cut -off = J s ("^)lcad.cut-off 



S 2 



2 f M^t L + i)-i , n^in^ + i) i 



= W(2^- + 4 ^ 7+2lh)' ^ 

These are the cross sections obtained by introducing a Debye screening mass by hand into 



the vacuum formula (see Eq. (23)). As already mentioned, there is no way that magnetic 
screening is incorporated in these formulas. However, one could argue that the contribution 
from the magnetic sector to the total and transport cross section is small. This is to be 
checked here. Since any significant change in the total and transport cross sections crucially 
influences the predictions for all the global quantities and the signatures of the quark-gluon 
plasma (see f|]), it is important to compare our medium modified cross section (both electric 
and magnetic screening taken into account) with that of the Debye screened cross sections 
obtained by using Eq. (p3|). In Fig. || we plot the total cross section obtained by using the 
medium modified propagator including electric and magnetic sreening (see Eq. (|25|)) and the 
total cross section obtained by using the in-vacuum propagator with a Debye screening mass 
as a cut-off (see Eq. Q2"7|)). In this figure, we use the values of s and a which correspond to 
a "realistic" situation at RHIC It can be seen that for the values of s and a considered 
here, the difference in the two cross sections increases for decreasing values of s. In Fig.H, we 
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plot the transport cross sections for the two cases described above (Eq.fl26|) and (p8|)). The 
transport cross section directly enters into the collision term (see Eq. |5|) and so determines 
the evolution of the quark-gluon plasma. Similar to the total cross section, it can be seen 
that the medium modified transport cross section is enhanced in comparison to the other 
for smaller values of s. Hence, it can be expected that all the global quantities and the 
equilibration time might change when the vacuum propagator is replaced by the medium 
modified one. However, this has to be checked by including all these features in a self 
consistent transport study. For example, in the above plots, for simplicity we have used 
a constant value of a = 0.3. In "realistic" situations at RHIC and LHC the a, s, total 
cross section, transport cross section, as well as all other quantities are time dependent and, 
therefore, have to be calculated via the distribution functions of the partons, which in turn 
are obtained by solving the transport equation self-consistently Q]. As the self-consistent 
transport study involves extensive additional numerical work, the results incorporating the 
medium modified transport cross section will be presented elsewhere. 

In summary, the medium modified propagator has been used to evaluate the partonic 
scattering cross section for the process gg — > gg in the QCD medium by performing an 
explicit sum over the physical gluon polarizations. A magnetic screening mass from a non- 
perturbative study was used to show that the medium modified cross section is finite, diver- 
gence free, and is uniquely determined. The medium modified cross sections yield different 
results than these obtained by artificially introducing a Debye screening mass by hand into 
the vacuum formulae. This implies that the medium modified scattering parton cross sec- 
tions must be incorporated properly into the transport equations in order to study the 
production and possible equilibration of the minijet plasma at RHIC and LHC. Any change 
of the cross sections due to in-medium effects will crucially change the equilibration times, 
time evolution of the energy densities, number densities, temperatures of all degrees of free- 
dom and hence, all predictions of signatures of the quark-gluon plasma, both at RHIC and 
LHC. 
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Fig. 1 In the upper curve the difference of the total cross section in 
small angle approximation to the same quantity obtained by using the 
t-channel matrix element Icr^^^ — C(t-channei) I is plotted as a function 
of y/s (a = 0.3). The lower curve is the difference obtained for the 
transport cross sections. 
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Fig. 2 The total in-medium cross section and the total cross section 
regularized by putting a cut-off by hand ertot is plotted as a function of 
y/s. (a = 0.3) The first is significantly enhanced in comparison to the 
second. 
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Fig. 3 The in-medium transport cross section and the transport cross 
section regularized by putting a cut-off by hand <7 tra ns is plotted as a 
function of y/s. (a = 0.3) The first is significantly enhanced in compari- 
son to the second. 
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